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INVESTIGATION AND CALCULATIONS OF 
BASIC PARAMETERS FOR THE APPLICATION OF THE 
LASER DOPPLER VELOCIMETER 

By James F. Meyers 
Langley Research Center 

SUMMARY 

This paper is concerned with a recent development in a nonprobe technique for 
velocity measurement, namely, the laser Doppler velocimeter (LDV). The technique 
makes use of the Doppler effect on scattered laser radiation from small particles in the 
gas flow. Two LDV techniques currently in use, those being the single-beam technique 
and the cross-beam technique, are described and compared. The three modes of opera- 
tion of the cross-beam technique are compared with respect to the output characteristics 
and particle density requirements in the flow. The dual-scatter — forward-scatter mode 
is shown to be the most efficient. The dual-scatter — back-scatter mode, although not as 
efficient, has advantages which make the mode desirable for some applications. 

A computer program involving the basic LDV theory is presented and is used to 
derive the results for the comparison of the three modes in the cross-beam system. 

These results included the output signal-to-noise ratio and required seeding density as a 
function of the input laser power and scattering angle. By using a modified version of 
the program, the required input laser power was determined as a function of molecular 
density for a fixed signal-to-noise ratio, based on only molecular scatter. 

A test example is presented to show how the basic theory and the computer program 
are used to determine the operating criteria required for LDV operation in a particular 
flow in a wind tunnel at Mach 6. 


INTRODUCTION 

The aerodynamist and gas dynamist are very interested in turbulent gas flows and, 
in particular, in determining various turbulent-flow parameters such as the velocity vec- 
tor components. 

In the past few years a nonprobe technique has been under development by several 
research groups (as reported in refs. 1, 2, 3, and 4) that will measure the instantaneous 
velocity at a point in the flow. The technique involves scattering laser light from gas 



molecules and small particles in the flow. Due to the Doppler effect on the scattered 
light, the velocity vector components can be determined and the flow can be characterized. 
This system is generally called the laser Doppler velocimeter (LDV). 

Since the advent of the laser, much work has been done on the study of scattered 
laser light from gas molecules and small particles suspended in the gas. This work has 
been mainly concerned with atmospheric study using the Light Detection and Ranging 
technique (ref. 5). Work has also been performed on the development of theory used to 
describe the operation of the LDV. The purpose of this study is to investigate broadly 
the relation between the scattering of light from moving particles and the apparatus neces- 
sary to extract velocity information from the scattered light. Presently, the LDV charac- 
teristics are calculated independent of one another and related experimentally. This 
method can do little to predict the parameter requirements for LDV operation in a par- 
ticular flow. By relating all LDV characteristics, LDV systems may be theoretically 
applied to any flow under consideration. A computer program presented in this paper 
contains the related theory necessary to calculate all LDV parameters. The resulting 
plots of a study presented in the paper relate signal-to-noise ratio, photomultiplier output 
current and power, required seeding density, and spatial resolution with the LDV scat- 
tering angle. The scattering angle and the photodetector frequency response determine 
the upper limit of particle velocity that can be measured. Thus, all LDV characteristics 
are related and considered. 

The program has been applied to two LDV configurations currently in use. These 
are the single-beam system in the reference-scatter mode and the cross-beam system 
in the reference-scatter, dual-scatter — forward-scatter, and dual-scatter — back-scatter 
modes. As an example of the use of the program in determining the required param- 
eters for LDV operation, a typical gas flow is analyzed for possible application of an 
LDV to survey the velocity field. 


SYMBOLS 


A CL 

A r 

B 


The symbols 

[acl] 

[ar] 

00 

00 


in brackets are the equivalent symbols used in the computer program, 
cross-sectional area of receiver optics, meters^ 
area of receiving lens, meters^ 
statistical factor due to dynode emission 
velocity of light, 3 X 10^ meters/second 


2 



D 

[»] 

Da 

[DA] 

D r 

[DR] 

E r cos co r t 


E s cos w s t 


E t 


e 

H 

F 

M 

FI 

[FI] 

f 


f D 


f L 


fo 


f P 


fs 


g p 

M 

h 

W 


i(t) [XIOUT] 


diameter of input aperture, meters 
Airy disk diameter at sample volume, meters 
diameter of receiving lens, meters 
reference electromagnetic wave, volts/meter 
scattered electromagnetic wave, volts/meter 

total electromagnetic wave, E r cos o> r t + E s cos o> s t, volts/meter 

electron charge, 1.6 x lO - ^ coulomb 

laser power density at sampling point, watts/meter^ 

input lens focal length, meters 

depolarization attributable to anisotropy of scatterer 

light frequency observed by stationary observer of laser beam 
scattering from moving particle, hertz 

laser frequency, hertz 

light frequency observed by stationary observer positioned along 
velocity vector of moving light source, hertz 

light frequency observed by moving particle from stationary laser 
source, hertz 

source light frequency, hertz 
gain of photomultiplier 

Planck's constant, 6.626 x 10”^ joule-second 
total photomultiplier output current, amperes 
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W‘> 


dc component of photomultiplier output current, amperes 


i if (t) 


i^O!pn,0j [XIMEE] Mie intensity function, 


il(Q!l,n,0), 

i M (« 1 ,n,e) 

<‘ lf 2 > 

<‘n 2 ) 

k 

L 

N 


N 


Mie 


N t 


ac component of photomultiplier output current, amperes 

M) + i ||( a 1 » n ,g) 


[XK] 

[xl] 

Oxn] 

[xnmie] 


intensity of light with electric vector perpendicular and parallel 
to plane through direction of propagation of incident and 
scattered beams, respectively 

mean-square photodetector response, amperes^ 


mean-square noise response in detector due to shot noise, 

9 

amperes^ 

2 7T — 1 

wave number of incident radiation, — , meter 

A 


SC 


length of sample volume, meters 

number of photons per second impinging on photocathode surface, 
photons/second 

number density of particles at sample volume, 
particles /meter ^ 

number density of molecules at sample volume, molecules/meter^ 

[sRE f) index of refraction 

[pq] input laser power, watts 

unscattered laser power used as reference signal, watts 

[XPOUT] output power of photomultiplier, dBm 

[PS] scattered laser power used as signal, watts 

[PSC] laser power scattered from flow in acceptance cone of 
receiver, watts 



q 


attenuation of wave (transmissivity) in propagating from 
scatterer to detector lens 

R 

[ R ] 

output impedance of matching unity gain amplifier on photo- 
multiplier, ohms 

r 


radius of particle, meters 

S l 

[slast] 

final particle diameter, microns 

Ss 

1 1 

starting particle diameter, microns 

S/N 


signal -to-noise ratio 

(S/N) dB 

M 

signal -to-noise ratio, decibels 

t 


time, seconds 

V s 


particle velocity, meters/second 

Z' 

[z] 

distance from receiving lens to sample volume, meters 

a 


angle between particle velocity vector and laser beam, degrees 

a l 

OO 

2 jr r 

particle size parameter, — 

a 


average polarizability of air, meter s^ 

1 3 


proportionality constant of photomultiplier 

A 


depolarization factor 

AS 

[delta] 

incremental particle diameter, microns 

AV 

[V] 

sample volume in flow, meters^ 

A6> 

[deth] 

incremental cross-beam angle, degrees 

Af 

[delf! 

noise bandwidth, hertz 
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Afi 


D 


Doppler frequency, hertz 


V 

e 


o' 

A 


CT Mie 


$ 


01 

0 2 

ft 


transmission coefficient of receiver optics 
jjELl] transmission coefficient of receiving lens 
[EL2j transmission coefficient of beam splitter 
[XNU] quantum efficiency of photomultiplier 

angle between unscattered laser beam and signal optics, degrees 
[dTHELA] final cross-beam angle, degrees 
[DTHETA] starting cross-beam angle, degrees 

angle between laser beam and scattered beam, degrees 
QxLAmJ wavelength of laser beam, meters 


optical frequency of laser, hertz 


Mie scattering cross section per unit volume, meter 


-1 


-1 


Rayleigh scattering cross section per unit volume, meter 

angle of rotation of receiver optics about sample volume in 
cross-beam system, degrees 

angle between velocity vector and beam 1 in cross-beam sys- 
tem, degrees 

angle between crossed laser beams, degrees 
solid angle viewed by signal optics, steradians 
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co r angular frequency of reference electromagnetic wave, 

radians/second 

w s angular frequency of signal electromagnetic wave, radians/second 

COMPARISON BETWEEN THE SINGLE -BEAM 

And cross-beam systems 

The laser Doppler velocimeter as developed by Foreman, George, and Lewis 
(ref. 1) has been used to measure accurately high-velocity flows in wind tunnels without 
disturbing the flow as with previous techniques (ref. 2). This system, which will be 
referred to as the reference-scatter single-beam system, employs a focused laser beam 
which is scattered from particles in the gas flow. A portion of the scattered light is 
sampled and mixed with a portion of the unscattered laser beam and collected by an 
optical-photomultiplier system where the two light beams heterodyne to yield the differ- 
ence frequency between the two light beams. This difference frequency, or Doppler fre- 
quency, is related to the particle velocity in the flow. Particles of micron size are 
expected to follow the gas velocity within a few percent and, for the purposes of this study, 
are assumed to do so. The relationship between the particle velocity and Doppler fre- 
quency, as developed in appendix A as equation (A10), is 

Afjj = -^•|cos(q' - 6') - cos aj 

where 

Afj} Doppler frequency, Hz 

V s particle velocity, m/s 

X wavelength of laser beam, m 

a angle between particle velocity vector and laser beam, deg 

d' angle between laser beam and scattered light, deg 

A typical reference -scatter single-beam LDV is shown in figure 1. 

A new LDV technique has recently been developed (ref. 3) which involves an optical 
system that divides the input laser beam into two parallel beams which are then focused 
to a point, thus crossing the two beams. If a particle is located at the crossover point, 
the particle will scatter light from both beams. By placing a photomultiplier along the 
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axis of one of the beams, a reference-scatter system is constructed, since the photo- 
multiplier also collects scattered radiation from the other beam. Since the scattered 
radiation from the other beam radiates in all directions, part of this scattered light is 
coaxial and parallel with the reference beam and thus automatically alines the photo- 
multiplier. This technique is referred to as the cross-beam system and is shown in 
figure 2. 


If the receiving optics for the cross-beam system were located at a position mid- 
way between the two beams, the optics would then collect scattered light from both laser 
beams. Since the scattered light collected from each beam is of a different frequency, a 
heterodyne frequency is obtained. This configuration is shown in figure 3. From the 
theoretical development in appendix A, it is shown that the heterodyne frequency which is 
equal to the Doppler frequency is independent of the placement of the receiver optics, as 
long as the receiver optics are in the plane of the beams. The Doppler frequency thus 
derived (eq. (A18)) is 


Af D = 



where 


Afj) Doppler frequency, Hz 

V s particle velocity, m/s 

X laser wavelength, m 

angle between two laser beams, deg 


It can now be seen that the dual-scatter technique has an advantage over the 
reference-scatter technique. By referring to equation (A10), it is found that the Doppler 
frequency is dependent upon the angle of viewing the scattered radiation O' . Thus the 
Doppler frequency will vary across the viewing lens since the finite lens views a band of 
scattering angles. This problem is not found with the dual-scatter technique. As shown 
in equation (A18), the Doppler frequency only depends on the angle between the input laser 
beams 0g anc * is thus independent of the finite acceptance angle of the lens. In order to 
reduce the instrument broadening error, the reference-scatter systems must use only a 
small portion of the collecting lens; however, the dual-scatter system, with no instrument 
broadening, may use as large a collecting lens as practical because, as seen from equa- 
tion (A20), the contribution from the lens element outside the plane of the cross beams is 
of second order. As shown later this difference is a major advantage of the dual-scatter 
technique. 
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Alinement problems developing from the complex receiver optical system in the 
single-beam technique constitute an undesirable trait. The cross-beam technique has 
simple optical systems that are easily alined which, from this standpoint, makes the 
cross-beam system the more desirable LDV technique. Also, as shown in figures 2 
and 3, the cross-beam system may be used in either the reference-scatter mode or the 
dual-scatter mode. Since the final analysis to determine the most efficient system is 
yet to be made, the single-beam system is not discussed further inasmuch as it is limited 
to only the reference-scatter mode and the characteristics of the reference-scatter mode 
are identical in the single-beam and cross-beam systems. 

SYSTEM THEORY 

Derivation of the Relation Between Signal-to-Noise Ratio 
and Scattered Laser Power 

The instrumentation used to process the photomultiplier output current containing 
the Doppler frequency has input signal-to-noise and input power requirements. Since 
the power requirements may be met by using amplifiers, the signal-to-noise requirement 
is the limiting factor. Thus, a relation between signal-to-noise ratio and collected 
scattered light power must be determined. 

If the reference-scatter mode is considered, there are two parallel and coaxial 
electromagnetic waves, the unscattered reference wave and the scattered signal wave, 
impinging upon the photocathode surface. Just before contact, the total wave may be 
described as the sum of the two waves as follows: 

Et = E s cos u) g t + E r cos o> r t (1) 

where 

Et total wave 

E s cos o) s t signal wave 

E r cos o> r t reference wave 

Since the photomultiplier is a square-law detector, the current i(t) at its terminals 
is given by the equation 

i(t) Ej. 2 = Eg^cos^OJgt + E r ^cos^o> r t + EgE r cos(u> s - o> r )t + E g E r cos(o> s + w r )t (2) 
with 100 percent heterodyne efficiency assumed. 
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Taking a time average of equation (2) , where the interval is short compared with 
the Doppler frequency but long compared with the optical frequencies, gives 


i(t) = (3 


E 2 E r 2 

+ Eg E r cos (w s - w r )t 


( 3 ) 


where j3 is a proportionality constant of the photomultiplier. The photomultiplier is 
incapable of generating currents at optical frequencies 2o; s , 2w r , and w s + w r> so 
these terms do not appear in the output. The output current is the sum of two components: 
the dc current i^ c (t) and the Doppler frequency current i^(t) expressed as 


i dc (t) = 0 




(4) 


i if (t) = /3E s E r cos(o> s - w r )t 


(5) 


Thus, the Doppler frequency may be determined from the output of the photomultiplier. 
To continue the derivation: 


i(t) = 


1 + 


i if( t )~ 

W*) 


l dc 


(t) 


or 


i(t) = 


2E g E r cos(cj s - w r )t 

1 "1* n n 


E s 2 + E r 2 


dc 


(t) 


( 6 ) 

(7) 


Thus 


i if (t) = 


2E g E r cos(a> s - w r )t 


Eo 2 + E r 2 




( 8 ) 


The mean-square photodetector response ^-if / > where the brackets represent a 
time average over an interval of 1 Doppler frequency cycle, can then be computed from 
equation (8) and can be expressed as follows: 


/. 2\ _ 2P s P r ^c 2 ® 
\ if /~ (P s + P r ) 2 


(9) 


The noise sources for this system are considered to be composed of the dark- 
current shot noise, signal shot noise, background shot noise, and thermal noise (Johnson 
noise). Since the background light level can be controlled by the experimenter to a large 
degree, shot noise derived from the background light level is not considered in this paper. 
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It can be shown that in normal LDV operation the signal shot noise will be at least two 
orders of magnitude greater than either dark-current noise or Johnson noise. There- 
fore, it is assumed that this is the case in the present paper and that the total noise factor 
is due to signal shot noise only. Thus, the mean-square noise response <^i n ^ in the 
detector is given by 

^n ^ = 2e AfBgp i(j c (t) (10) 

where 

e electron charge 

Af noise bandwidth 

B statistical factor due to dynode emission of shot noise 

g D gain of photomultiplier 


A signal -to-noise ratio is defined as 

S = <Hl) = 2P s P r i dc 2 (t) 

N ^n 2 ) 2 ( p s 2 + 2p s p r + p r 2 ) e AfB Sp W 4 ) 

or 


S_ _ p s p r McW 

N_ ( P s 2 + 2P s p r + P r 2 ) e gp AfB 

but 


(ID 


( 12 ) 


6g P 


idc (t > 

( p s + p r)^ 


(assuming complete photoelectron collection). Thus 

S = ^sPr^s + Pr) _ ?7 p s p r 

N (P s 2 + 2P g P r + P r 2 ) AfBhf ( p s + P r ) AfBhi/ 


(13) 
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where 


77 quantum efficiency of photomultiplier 

h Planck's constant 

v optical frequency of laser light 


If P r is kept much larger than P s , equation (13) reduces to 
S _ 

N fm AfB 

or, with S/N given in dB, 

(2) = 10 log — 

\N/dB h v AfB 


(14a) 


(14b) 


In order to determine the ratio P r /P s when this condition applies, equation (13) is 
plotted with S/N as a function of P r for a constant P s . The constants in equation (13) 
for this test are based on a laser wavelength of 5145 A, an S-20 photomultiplier response, 
and a 1.0-MHz noise bandwidth. The resulting plots are given in figure 4. It may be seen 
that a P r /P s °f 4 gives a signal-to-noise ratio within 1 dB of the maximum signal-to- 
noise ratio. Therefore, if the ratio P r/ /P s is kept at 4 or greater, then equation (13) 
will reduce to equation (14a). 


If the dual-scatter system is considered with the optics located midway between the 
two laser beams, either forward-scatter or back-scatter, the signal wave and the refer- 
ence wave will both be scattered waves and of equal power. Since P s = P r in 
equation (13) 


S _ ^s 2 

N 2P S AfBhy 


(15) 


or 

S _ V * s 

N 2 AfBlm 


(16a) 


or, with S/N given in dB, 


^ = 10 log — 

,N/dB 2 AfBlm 


(16b) 
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Thus the relationships between signal -to-noise ratio and scattered power have been 
developed for both the reference-scatter system and the dual-scatter system. 

Rayleigh and Mie Scattering Functions 

From classical electromagnetic theory, the scattering process is a function of the 
particle size relative to the wavelength of the impinging light and the optical characteris- 
tics of the material of the particles. In the case of scatter from molecules 2r « X the 
theory has been developed by Rayleigh in reference 6. For particles on the order of the 
wavelength in size, the appropriate theory is given by Mie in reference 7. 

If a laser beam is focused to a small sample volume, the scattered power received 
by the detector lens is given by the equation (ref. 5): 

P sc = Fq(<T r + cr M i e )aTO (17) 

where (a r + is the macroscopic scattering cross section per unit volume for 

Rayleigh scatter and Mie scatter, respectively. The factor q is the attenuation of the 
wave (transmissivity) in propagating from the scatter er to the detector lens. The quantity 
AV is the sample volume as determined by the volume common to the crossing laser 
beams. The factor is the solid angle subtended by the receiving lens, and F is the 
laser power per unit area at the sample volume. 

The macroscopic Rayleigh cross section for forward or back scatter is, as obtained 
from reference 6, 

u r = t 1 + co - s . 2 -g ) k 4 a 2 N r f (18) 

2 


where 

k 

a 

N r 

f 

A 

e 


2tt 

wave number of the incident radiation, — 

A 

polarizability 
molecular number density 

depolarization attributable to the anisotropy of the scatter er, ^ 

6 - 7 A 

depolarization factor 

angle of scattered light received by the viewing optics 
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Based on the assumption that the particulate matter present in the flow may be con- 
sidered a collection of spherical particles of index of refraction n, the absolute Mie 
scattering cross section is given by the equation, as obtained from reference 7: 


CT Mie " k 2 Z ^"l’H’^Mie 

where 

+ i|i(ai,n,e) 

i^Q!j,n,0J Mie intensity function, — m L 


il^pM), 

i, 1(0-1, n,e) 


intensity of light with electric vector perpendicular and 
parallel to the plane through the direction of propagation 
of the incident and scattered beams, respectively. 

(These intensity functions were calculated by using the 
computer program developed by McCormick in 
reference 5.) 


a l 


particle size parameter, 


2rr 

X 


N Mie 


number density of particles at the sample volume 


(19) 


From equations (17), (18), and (19) the molecular and particulate number densities 
are related with the input laser power to the required scatter laser power. 


Expansion of Basic Theory 

In presenting the theory to obtain the criteria for the application of the laser 
Doppler velocimeter, consideration must be given to the limitations placed on the system 
by the input requirements for the processing instrumentation that will be used to condi- 
tion the Doppler frequency signal. These requirements are a minimum signal-to-noise 
ratio and an input power range. Since the power may be adjusted by means of amplifiers 
or attenuators, the signal-to-noise ratio becomes the major requirement. 

From equation (16b), the scattered laser power is related to the signal-to-noise 
ratio for the dual-scatter LDV. Solving for the scattered or signal power yields 

Ps = 2 AfBhl^lO^ (20) 

Taking the optical losses in the receiver into consideration gives 
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( 21 ) 



where e is the transmission coefficient of the receiver optics. 


Equation (17) then relates the laser power scattered from the flow to the scattering 
functions and thus to the molecular and particle densities. This equation, repeated for 
the convenience of the reader, is 


Psc = F< 1 ( CT r + CT M ie) AV ^ 

where F is the laser power density at the sample volume and is defined by 


F = 


qP 

A CL 


( 22 ) 


where q is the transmissivity of the medium and P is the input laser power, 
cross-sectional area of the sample volume A^l is defined 


a CL “ 


^ D a 2 

4 


The 


(23) 


where D a is the diameter of the sample volume and may be approximated by the Airy 
disk diameter (ref. 8), which is expressed as 


^ 1.22X FI , 

Ua D { 

where X is the wavelength of light, FI is the input lens focal length, and D is the 
input aperture diameter. 

The sample volume AV is in the form of an ellipsoid, shown in the following 
sketch, as described by the Airy disk diameter of the focused laser beam: 
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This volume may be determined by the Airy disk diameter of the focused laser beam as 
based on the input focusing lens and the sine of the angle between the two laser beams as 


AV = 


* p a 3 
6 sin <p 2 


(25) 


It is assumed that the entire sample volume is viewed by the receiving optics. If a larger 
volume is viewed by the receiving optics only the dc level will increase, with no change 
in spatial resolution with respect to the Doppler sample volume. This is because only 
scattered light from the beams scattering simultaneously from the same particle can com- 
bine in order to obtain the heterodyne effect. 

The solid angle viewed by the receiving optics is given by 
A r 

ft = ~ (26) 

Z 2 

where A r is the area of the receiving lens and Z is the distance between the sample 
volume and the receiving lens. Thus the equations for determining the necessary criteria 
for application of a LDV system are known. 


APPLICATION OF BASIC THEORY 
Development of a Computer Program 

By programing the basic theory for the laser Doppler velocimeter, a design study 
may be made over a variety of physical parameters. Also, effects of instrument varia- 
tions may be determined. 

Consider a program that is to yield the necessary criteria for LDV operation in 
low-density wind tunnels, for example, a Mach 6 wind tunnel with a stagnation pressure of 
4.14 MN/m^ (600 psi) and a stagnation temperature of 204° C (400° F). With a closer 
look at the basic theory, it is found that for these conditions, if particles must be added 
to the flow in order to obtain the necessary scattered laser power, the molecular scat- 
tering contribution may be neglected. From appendix B, considering the example flow, 
the Rayleigh scattering cross section cr r is 2.27 X 10"^ m~^. If 1 particle /cm^, with a 
particle diameter on the order of a micron, is added to the flow, the Mie scattering cross 
section is 1.7 x 10“^ m”^. For this condition the molecular scattering contribu- 

tion may be neglected. For flows where the molecular densities are higher, this may not 
be the case. The molecular scattering results presented later give some approximation 
of the regions where molecular scattering must be considered. 

The computer program containing the basic theory is presented in appendix C. In 
this program the contribution of scattered power due to molecular scatter has been 
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neglected. The program as presented in appendix C is for the problem in which the 
required signal-to-noise ratio is an input and the minimum particle density necessary to 
obtain the required signal-to-noise ratio is determined. The LDV used in this case is 
the cross-beam dual-scatter system. The program will accept either forward-scatter or 
back-scatter problems. With a few modifications the program may be converted to handle 
a reference-scatter system. The program may also be converted for solving problems 
with particle density as an input and resultant signal-to-noise ratio determined. 

Results of Application of the Program 

The cross-beam system is the LDV considered. The system is used in the three 
modes of operation so that a comparison may be made between the modes. The modes 
are described as follows: 

(1) The dual-scatter — forward-scatter mode with the receiver lens located between 
the cross beams. The receiver lens is moved as the cross-beam angle is varied in such 
a way that the lens just remains between the beams. 

(2) The dual-scatter — back-scatter mode with the input focusing lens used as the 
receiving lens 

(3) The reference-scatter mode with the receiver located so that the distance 
between the sample volume and the receiving lens is the same as the lens distance for 
the dual-scatter — forward-scatter system with a cross-beam angle of 4° 

Each mode has its advantages as described previously. By the use of the program, a 
quantitive comparison is made of the modes as the cross-beam angle is varied from 1° 
to 8°. The cross beams for the dual-scatter modes will be of equal power whereas the 
cross beams for the reference-scatter mode will be 99 percent and 1 percent with the 
lower power beam used as the reference beam. The power ratio in the reference -scatter 
mode is used to obtain the reference-to-signal power ratio needed to maximize the signal- 
to-noise ratio as described previously. Also, by keeping the reference-beam power low, 
the photomultiplier is kept out of saturation. 

In order to perform the calculations in the program, a few basic criteria must be 
fixed. For an example, the chosen criteria are 

Laser: 1.0-W Argon at 5145 A (unless specified otherwise), 1.5-mm beam diameter 

Input focusing lens focal length: 0.6096 m (24 in.) 

Receiving lens diameter: 0.0762 m (3 in.) 

Photomultiplier response: S-20 

Photomultiplier gain: 8.0 X 10^ 

Noise bandwidth: 10® Hz 
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The seeding particles are assumed to be uniformly distributed, having diameters 
from 0.4 /z to 2.0 /a, with an index of refraction of 1.48. 

In figure 5, the sample volume is determined by the input focusing lens, which is 
the same for all three modes, as based on the Airy disk diameter. As is shown, the vol- 
ume decreases as the angle between the cross beams increases. The length of the sample 
volume is given in figure 6. The sample-volume length also decreases with an increase 
in the angle between the cross beams. Thus the spatial resolution of the LDV is deter- 
mined as a function of cross-beam angle. 

For the fixed inputs given, the signal -to-noise ratios, with a seeding particle density 
of 1000 particles/cm^ and for cross-beam angles from 1° to 8°, have been computed. The 
plots for the three modes of operation are given in figure 7. Taking into consideration 
the placement of the receiver optics and the theoretical discussion, a closer examination 
of these plots may be made. 

For a large signal-to-noise ratio, the dual-scatter — forward-scatter mode is found 
to be the best of the three modes. The primary reasons for the large signal-to-noise 
ratios and the shape of the curves may be found by referring to equations (A18), (16b), 
and (17) and figure 5. Equation (A18) states that the dual-scatter mode contains no instru- 
ment broadening; thus, a large receiving lens may be used, which allows a large viewing 
solid angle. By keeping the receiving lens just between the two beams, the solid angle 
viewed increases with increased cross-beam angle, which allows increased collection of 
scattered power. On the other hand, the sample volume decreases with an increase in 
cross-beam angle and thus limits the number of scattering centers. Also, the signal-to- 
noise ratio is less by a factor of two for the dual-scatter modes than for the reference- 
scatter mode. By taking these parameters into consideration in equation (17), the scat- 
tered power received by the receiving optics is determined. It can be seen in figure 7 that 
the solid angle increasing with increased cross-beam angle more than compensates for 
the loss of signal-to-noise ratio in equation (16b) and the decrease in sample volume. 

For example, the sample volume decreases by a factor of five with an increase in cross- 
beam angle from 1° to 5°, whereas the viewing solid angle increases by a factor of 25. 

The dual-scatter — back-scatter mode has the same advantages and disadvantages 
as the dual-scatter— forward-scatter mode except that the viewing solid angle is fixed 
since the input lens is used as the receiving lens. This is the reason the curve decreases 
with increased cross-beam angle since the decreasing sample volume is the dominating 
factor in equation (17). Also, the Mie cross section is a factor of about 10 less in the 
back-scatter mode than in the forward-scatter mode, which reduces the signal-to-noise 
ratio and output power. 

The low effective signal-to-noise ratio obtained from the reference -scatter mode is 
due to the fact that the useful viewing solid angle is determined by the reference -beam 
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cross section at the receiving lens. Thus, the scattered power as determined by use of 
equation (17) is very low and decreases with the sample volume decrease for increased 
cross-beam angle. The scattered power is so low that the compensation due to a more 
efficient signal -to -noise relation does not bring the effective signal-to-noise ratio up to 
the levels of the dual-scatter modes. 

Figure 8 is obtained by performing calculations to determine photomultiplier output 
current based on LDV operation plotted in figure 7. These calculations are made by con- 
verting the input laser power to the photomultiplier from watts to photons/second and 
applying the following equation: 

i(t) = Ng p er? (27) 

where 

i(t) photomultiplier output current, A 

N number of photons per second impinging on the photocathode surface 

gp gain of the photomultiplier 

e electron charge, C 

7 } quantum efficiency of the photomultiplier for the laser wavelength 

With the output current known, the output power may be determined for a 50-ohm 
characteristic impedance coaxial transmission line. With an output power in the form of 
dBm, the results for the three modes are compared in figure 9. 

By working backwards with the signal-to-noise ratio as the input limit, the minimum 
particle density that will satisfy the instrument requirements (for example, 15 dB) may 
be calculated. The particle densities thus obtained are given in figure 10. On the basis 
of required minimum seeding particle density, the dual-scatter — forward-scatter mode 
is again the most favorable. In fact, it would seem fromthe curve that at a slightly 
larger cross-beam angle there would be enough scattered power from 1 particle/cm^ to 
obtain an acceptable signal-to-noise ratio. However, the signal obtained would be very 
intermittent and is not discussed further since a continuous output is assumed in this 
paper. 

In figure 11a comparison of the dual-scatter — forward-scatter mode and the 
reference-scatter mode is made for two different lasers, a 50-mW He/Ne and a 1.0- W 

o 

Argon at 5145 A. It was found quite dramatically just how much more efficient the dual- 
scatter technique is than the reference-scatter technique since for a cross-beam angle 
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above 2.5° the required particle density for the use of the dual-scatter mode with the 
He/Ne laser is less, even up to a factor of 10, than the reference -scatter mode with the 
Argon laser. 

For applications where it is impossible or undesirable to add seeding particles, it 
is worthwhile to consider the possibility of Rayleigh scattering from the molecules them- 
selves. The computer program is now modified to handle molecular scattering param- 
eters. The portion of the program containing the calculations for determining Mie scat- 
tering parameters is removed and replaced with the calculations for Rayleigh scatter. 

By using the most efficient LDV (the dual-scatter — forward-scatter system) with a cross- 
beam angle of 6°, the input laser power requirements may be found for particular molec- 
ular densities (based on the 5145 A Argon line) that would satisfy the 15-dB signal-to- 
noise ratio in the example. The plots of molecular density as a function of input laser 
power are presented in figure 12 for three gases: helium, air, and carbon dioxide. 

To illustrate laser wavelength dependence on Rayleigh scatter, consider the same 
test with air and vary the laser wavelength from 3000 A to 7000 A in steps of 500 A. The 
results are given in figure 13 for an S-20 photomultiplier response. The most efficient 
response is around 3500 A. In order to determine the possibilities of using lasers cur- 
rently available, figure 14 is presented to illustrate the same computation with specific 
laser wavelengths. The lasers used are 


Argon (frequency doubled) 

2573 A 

ultraviolet 

Helium-Cadmium 

3250 A 

ultraviolet 


4416 A 

blue 

Argon 

4880 A 

blue-green 


5145 A 

green 

Helium-Neon 

6328 A 

red 

Ruby 

6943 A 

red 


The CC >2 laser (wavelength, 10.6 pm) is in the infrared and is not detectable with the S-20 
photomultiplier response. 

It may be seen in figure 14 that the Helium-Cadmium laser is the most efficient, 
but at present the maximum power obtainable is about 50 mW, which is too low to be con- 
sidered. The next most efficient is the Argon laser which does have the required output 
power in experimental lasers that would satisfy the requirements for Rayleigh scattering 
of air with 1 atmosphere pressure at 23° C. The ruby laser also has the required output 
power even with its long wavelength, but it is a pulsed laser which does not allow con- 
tinuous measurement of velocity. 
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Thus with a high-power Argon laser, the velocity of a gas flow of 1 atmosphere 
static pressure may be measured directly. If this is done, Brownian motion of the mole- 
cules must be taken into consideration when reducing the resulting data. 

Determination of Operating Criteria for an LDV in a Test Flow 

Consider a flow with a mean core velocity of 1200 m/s and a core diameter of 
15 cm, with a 7.5-cm-thick boundary layer about the core. It is desirable to scan the 
velocity field in this flow. Thus the mean velocity range is 0 to 1200 m/s. A turbulence 
study is also desirable so that the final velocity scan will contain a turbulent velocity 
deviation about the mean velocity. In transposing the velocity measurements to Doppler 
frequency there will be a mean Doppler frequency representing the mean velocity, a 
Doppler frequency deviation representing the magnitude of the velocity deviation, and a 
rate or modulation of the Doppler frequency representing the rate at which the velocity 
changes. Thus the output Doppler frequency is analogous to an FM radio signal where 
the mean frequency is the FM carrier frequency, the deviation is the amplitude of the sig- 
nal, and the modulation is the frequency of the signal carried by the carrier frequency. 

The tunnel windows in this example are taken to be 4.0 cm in diameter and to be 
located along an axis perpendicular to the flow. 

To begin the LDV design, the LDV restrictions previously used in the computer 
program will be adequate. The conditioning instrumentation used will be similar to the 
units used at the George C. Marshall Space Flight Center as discussed in reference 2. 
This instrumentation is based on a wide -band frequency tracker with automatic tracking 
and compression of signals. This tracker is considered in this example because it is 
typical of the readout instrumentation used with an LDV. The spectrum analyzer readout 
is not considered since it does not yield information that is readily reduced. 

The frequency tracker requirements are 

Input power: 0 dBm to -50 dBm 

Input signal -to -noise ratio: greater than 15 dB (1-MHz noise bandwidth) 

Frequency response: 5 MHz to 200 MHz 

Frequency deviation: 0 to 50 percent of carrier frequency or 30 MHz, 

whichever is lower 

Frequency modulation: 0 Hz to 100 KHz 
Thus the LDV restrictions are complete. 

Because of the window restriction of a 4.0-cm diameter, the greatest cross-beam 
angle that may be used is 7.5° for a flow diameter of 30 cm. The maximum cross-beam 
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angle as restricted by the frequency response of the conditioning instrumentation may be 
calculated from equation (A18) which is 


Af D 



and can be written 


-1 Af D A 

’2 = 2sin -2VT 


The mean velocity, 1200 m/s, may be measured with a 4.9° cross-beam angle. In order 
to allow for velocity fluctuation, the cross-beam angle is reduced to 4°. 

With these two limitations the cross -beam-angle range is fixed. From the fre- 
quency response limitations of the conditioning instrumentation, the velocity ranges for 
the cross -beam-angle limits were determined from equation (A18) and are shown in the 
following table. The deviations shown in the table were determined by the conditioning 
instrumentation. 


Cross-beam 
angle, (p 2 , deg 

Minimum velocity, 
m/s 

Maximum velocity, 
m/s 

Mean 

Deviation 

Mean 

Deviation 

4 

59 

±20 

1254 

±220 

7.5 

30 

±10 

668 

±118 


Since the 4° cross-beam angle gives the largest velocity range, the remainder of the 
criteria are based on this angle. 

From figure 5 it is found that the sample volume for this example is 1.25 x lO - *® m^. 
The length of the sample volume found in figure 6 is 5.2 mm. Thus the spatial resolution 
is determined. If a 1.0- W Argon laser is used, the required particle density is found 
from figure 10 to be 1500 particles/cm 3 for the reference-scatter mode, 130 particles/cm 3 
for the dual-scatter — back-scatter mode, and about 7 particles/cm 3 for the dual-scatter — 
forward-scatter mode. 

In order to measure the velocity cross section of the flow, the LDV must be trans- 
lated so that the sample volume, viewed by both the input and receiver lenses, moves 
across the flow. Since the dual-scatter — back-scatter mode has all components located 
on the same side of the tunnel, the system may be translated easier than the other modes. 
From all geometric and flow restrictions, the design engineer now knows the best LDV to 
use and the seeding requirements for that system. 
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Errors 


When the LDV requirements are met and a system is installed in a test condition, 
several sources or error must be considered. Also, requirements must be altered 
slightly to satisfy noise sources neglected in this paper. 

The accuracy of the final measurement is directly dependent upon the precision of 
the alinement of the LDV. In this consideration, errors may develop in the measurement 
of the angle between the cross beams, measurement of the laser-beam diameter entering 
the focusing lens, focal length of the lens, and location of the receiving lens. Although 
these errors are minor, a greater problem is found in that the spatial resolution is finite 
which can create large sources of error if the velocity gradient at the sample point is 
large. For example, with a 4° cross-beam angle, the sample volume is 5.2 mm long 
which means that there is a 3.4-percent error about a center velocity of 760 m/s for a 
velocity gradient of 50 m/s/cm. Another source of error is the accuracy of the condi- 
tioning instrumentation used to process the Doppler frequency; as an example, the fre- 
quency tracker has a 3-percent inaccuracy of measurement. 

The noise contributions that may be found are shot noise due to background light 
and noise introduced in the signal processing instrumentation. The net effect of including 
these sources of noise is to increase the level of the minimum detectable signal size. 

CONCLUDING REMARKS 

Two laser Doppler velocimeter (LDV) techniques, the single -beam and the cross- 
beam, have been described. The cross-beam system was shown to be the more versatile 
because it can be operated in three modes, the reference-scatter, the dual-scatter — 
forward-scatter, and the dual-scatter — back-scatter. Also, due to the simplicity of the 
cross-beam system it is easier to aline than the single-beam system. 

The basic theory for the calculation of the criteria for LDV operation was presented. 
A computer program incorporating the basic theory was also presented which allows the 
operating criteria to be determined easily for any LDV configuration described in this 
paper. The only parameters needed to make the calculations are tunnel and signal con- 
ditioning instrumentation requirements. 

By presenting an example set of restricting parameters, a series of operating 
criteria were determined. These results show that the dual-scatter — forward-scatter 
mode of the cross-beam LDV is the most efficient of the three modes. Although not as 
efficient, the dual-scatter — back-scatter system has two advantages: (1) ease of trans- 

lation since all components are located on the same side of the flow and (2) only one 
entrance window is required. 
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A typical design problem utilizing the techniques given in the paper is presented. 
The results obtained are the criteria necessary to obtain velocity measurements with the 
LDV. Thus, the design engineer may use the program to determine the optimum criteria 
for velocity measurement in a flow under study. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 10, 1971. 
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APPENDIX A 


DERIVATION OF DOPPLER FREQUENCY 


Consider a small particle moving along a velocity vector V s . If the particle 
radiates light energy and is viewed by a stationary observer located at a point along the 
velocity vector in front of the particle, the distance traveled by the light waves toward 
the observer in time t is (c - V s )t. The number of light wavelengths emitted in this 
time is f g t, where f s is the frequency of the source. Therefore, the frequency seen 
by the observer f Q is expressed (ref. 9) 


cf s t _ cf s 
0 = (c - Vs)t = c - V s 


(Al) 


However, if the light source is stationary and the observer is a particle moving away from 
the source, the frequency seen by the observer fp is 


f 


P 


= f L 



(A2) 


If the source is not along the velocity vector, the frequency seen by the particle is 
given by 


„ „ c - V q cos a 

h - *l h 


(A3) 


where a is the angle between the light beam and the velocity vector. 

As light is scattered from the moving particle, it now becomes the source. With 
the observer located as shown in the following sketch, the observer will see the fre- 
quency f D . 
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The frequency fp can be written 

£ _ £ £ 

D c - V s cos(Q! + Of P 

Substituting for fp in equation (A4) gives 

c - V c cos a c 


f D = 


:rf T 


c C - Vg COS(Q! + 6 ) ^ 


or 


c - V„ cos a 

f ® f 

^ c - V g cos(o! + 9) ^ 


(A4) 


(A 5) 


(A6) 


Reference-Scatter LDV 

If a reference-scatter system is used to determine the Doppler frequency, the angle 
relation is shown in the following sketch: 



Note that cos(o! + 6) in equation (A6) becomes cos (a + O') where O' is the angle 
between the observer optics and the reference beam. In the detector the two beams, the 
signal beam and the unscattered laser beam, are placed so that they are coaxial and 
parallel and impinge upon a photomultiplier where they heterodyne to yield the difference 
frequency between the light beams or the Doppler frequency. 

In order to obtain the Doppler frequency Afp, the source frequency f^ is sub- 
tracted from both sides of equation (A6) to give 


c - V Q cos a 

Af-pv = fj-v - f j = — f T - f -r 

D D L c _ Vg ^( 0 , _ e’) L L 


(A7) 
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or 


Af D - 


c - V s cos a - c + V s cos(a - O') 
c - V s cos(a - O') 


(A8) 


But c » V s cos(o< - O') for all velocities under consideration. Thus 


A e V s c cos^ - 0 T ) - cos a 

at D= X 

(A9) 

the Doppler equation becomes 


Afj) = — [coslo; - O') - cos aj 

(A10) 


Dual-Scatter LDV 

Now consider the dual-scatter system. If the observer is located as shown in the 
following sketch, the observer will collect scattered light from both beams. 


Beam 1 Beam 2 



If equation (A6) is applied to beam 1, the light frequency observed becomes 

f - c ~ v s cos a f 
^1 c - V s cos (a + 0) ^ 

where a and 0 are redefined by the geometry as 

Q!=0j O' + 0 = <f* 


(All) 
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Therefore 

c - V„ cos 0, 

f n = 2 i f T 

^1 c - Vg cos $ L 

If equation (A6) is applied to beam 2 , with a and 9 redefined as 

O' = 0j + 0 g G! + 0 = <f> 

the light frequency observed becomes 


(A12) 




c - V g cos (4^ + 0 2 ) 
c - V s cos $ 


f T 


(A13) 


The Doppler, or difference, frequency may then be found for the dual-scatter sys- 
tem by subtracting equation (A13) from equation (A12) to give 


c - V c cos 0 1 

Af-pv = fpv . - fpx — — - ~ 

D Dj D2 c - Vc cos <f> 


f L- 


c - V g cos (0 X + 0 2 ) 


c - V Q cos $ 


or 


Af n = 

u c - Vo cos 


— |c - V s cos - c + V s cos + 0 2 Jj 


(A 14) 


(A 15) 


But c » Vg cos for all velocities under consideration. Therefore 

*lVbi 


if D= c 


-|cos^0j + 0^ - COS 0^J 


(A16) 


Thus, the Doppler frequency is independent of the viewing angle of the observer for 
a dual-scatter system. This is not true for the reference -scatter system as shown in 
equation (A10). 

Reducing the Doppler frequency to a more usable form yields 


V 

Af D = -^-Jcos^j + 0 2 ) - COS j 

o ^2 

For the case of symmetry, 0^ = 90 — — 


Af V S 

Af D = —r— 


0 


cos (^90° — — + 0 2 ) - cos (90° — p 


0, 


Af 


D 


™s cos ( 90O J_2 


Af 2V s • 0 2 

Af D = -— sin-g- 


(A17) 


(A 18) 
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Thus the Doppler frequency conversion factor 


V s 



(A19) 


is only dependent upon the laser wavelength and the angle between the two input laser 
beams. Since it is desirable for the entire receiving lens to be used, the effect of col- 
lecting light out of the plane of the beams must be considered. The effect is found to be 
an additional cosine term in the denominator in equation (A15) so that the equation 
becomes 

Af D * c - v s cos# cos x [ c - V s cos *1 • c + V s cos (*1 + “>2)] < A20 > 


where X is the angle of scatter collected out of the plane of the laser beams. 

Thus a similar assumption that c » V s cos $ cos X is still valid, which means 
that the Doppler frequency as described by equation (A19) is valid over the entire col- 
lecting lens. 
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DETERMINATION OF THE EFFECTS OF RAYLEIGH SCATTER 

In order to determine the contribution of Rayleigh scatter and Mie scatter to the 
total scattered laser power, these parameters are calculated for a typical wind tunnel. 
The wind tunnel used in this example is a Mach 6 wind tunnel with a stagnation pressure 
of 4.14 MN/m2 (600 psi) and a stagnation temperature of 204° C (400° F). Based on the 
assumption that the flow is free of particulate material, only Rayleigh scatter from the 
molecules is considered. 

Begin by determining the Rayleigh scattering cross section from equation (18), 
which is 

(l -f cos 2 fl) k 4-2 N f 
2 

where 

k wave number, ^ or 1.22 x 10^ 

A 

a polarizability of air, 1.79 X 10" 30 m 3 

N r molecular density, 3.04 x 10 3 ^ molecules/m 3 

f depolarization factor (ref. 5), 1.054 

9 viewing angle of scattered light, 4° 

The result is a r = 2.27 x 10 - ^ m - *. 

In order to determine the effect of added particles, seed the flow with particles in 
a uniform distribution from 0.4 ju to 2.0 p in diameter so that the particle concentration 
is 1 particle/cm 3 . The Mie scattering cross section is found from equation (19) which is 

°Mie = “2 Z i (“l ,n,0 ) N Mie 
The result is cr Mie = 1.7 x 10“^ m“*. 

When the Mie scattering cross section is compared to the Rayleigh scattering 
cross section it is seen that the Rayleigh scattering cross section can be neglected if 
the particle concentration is 1 particle /cm 3 or greater if the molecular density is 
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3.04 X 10^ molecules/m^ or less. If higher molecular densities are in the flow, this 
may not be the case and these calculations must be performed to determine the contribu- 
tion of each scattering cross section. 
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COMPUTER PROGRAM FOR CALCULATION OF LDV PARAMETERS 


The organization of the computer program used to obtain the required parameters 
for a given LDV configuration is illustrated in the following flow chart: 
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A detailed explanation of the numbered blocks in the flow chart is as follows: 
Block 1: Input parameters needed 

(S/N)^ signal -to-noise ratio, dB 

P input laser power, W 

FI input lens focal length, m 

D r diameter of receiving lens, m 

Z distance from sample volume to receiving lens, m 

X wavelength of laser beam, m 

D diameter of input aperture, m 

0 S starting cross-beam angle, deg 

A 8 incremental cross-beam angle, deg 

81 final cross-beam angle, deg 

S s starting particle diameter, microns 

AS incremental particle diameter, microns 

final particle diameter, microns 
n index of refraction of particles 

Block 2: Printed output of input parameters 

Block 3: Calculation of the required scattered power impinging on the photomultiplie 

to satisfy the input signal -to-noise parameter from equation (2) 

(s/m 

_ 2hc AfBICr 10 ' 

S - A 77 
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Block 4: 


Block 5: 


Block 6: 


Block 7: 


Block 8: 


where 

h 


Af 

B 

V 


Planck’s constant, 6.626 x 10“34 j_ s 
velocity of light, 3 x 10® m/s 
noise bandwidth, 10 Hz 

statistical factor due to dynode emission of shot noise, 1.2 

quantum efficiency of the photomultiplier, 0.13 for 5145 A 
laser wavelength and S-20 response 


Calculation of the area of the receiving lens 
2 


Ay = 


7TDi 


(Cl) 


Calculation of the Airy disk diameter as based on the input focusing lens by 
use of equation (24) 




1.22A FI 
D 


Subroutine to calculate the Mie intensity function for a uniform distribution of 
particle size as given in the input parameters. 

Calculation of the wave number 

k = ^ (C2) 

Calculation of scattered power required from the flow (eq. (17)) 



where 

percent transmission of light through receiving lens, 0.92 
eg percent transmission of light through beam splitter, 0.45 
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Block 9: 


Block 10: 


Block 11: 


Block 12: 


Block 13: 


Block 14: 


Block 15: 


Block 16: 


Calculation of viewing solid angle by use of equation (26) 
A r 


a = 


Z 2 


Calculation of cross-sectional area of the sample volume by use of 
equation (23) 
ifD a 2 

a CL _ 4 

Calculation of sample volume by use of equation (25) 

3 


AV = 


7tD ? 


6 sin e 

Calculation of the length of the sample volume 
L = j 2Da 

sin $ 

Calculation of laser power density at sample volume 
P 


F = 


2A 


(C3) 


(C4) 


CL 


The factor of 2 is included because of the division of the input laser beam 
into two crossing beams for the dual-scatter system 

Calculation of number of photons per second impinging upon the photocathode 
surface 

2AP C 


N = 


he 


(C5) 


Calculation of the required number density to satisfy the input parameters 


N Mie " 


k 2 P 


sc 


F AVO 


I ‘(“l-"’ 9 ) 


(C6) 


Calculation of output current from the photomultiplier by use of equation (27) 


i(t) = Ng p er7 
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Block 17: 


Block 18: 


Block 19: 
Block 20: 


Printed output of calculated parameters 

P g required scatter power 

P gc scattered power required from flow 

D a Airy disk diameter 

Acl cross-sectional area of sample volume 

F laser power density at sample volume 

^ i^Ofpn,#) Mie intensity function 

r 

L length of sample volume 

AV sample volume 

ft viewing solid angle 

^Mie required particle density 

Calculation of output power from the photomultiplier 

P Q = 1000 i 2 (t) R (C7) 

where 

R output impedance of matching unity gain amplifier, 50 ohms 

(The 1000 is a conversion from watts to milliwatts.) 

Conversion of output power to output power in dBm 

Printed output of calculated parameters 

N number of photons per second impinging on the photocathode 

surface 
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Block 21: 
Block 22: 
Block 23: 


i(t) output current from photomultiplier 

P Q output power from photomultiplier 

IF statement to determine whether scattering-angle scan is complete 
Scattering-angle scan is complete; STOP 

Scattering-angle scan is not complete; increment angle and go to input 


A printout of the program is given on the following pages. 
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PROGRAM LASER( I NPUT , OUTPUT , T A PE 5= I NPUT , T APE6 = 0UTPUT , T APE 9 ) 
DIMENSION SREF ( 2 ) 


C THE PROGRAM WILL CALCULATE THE REQUIRED PARTICLE DENSITY* AND 

C OUTPUT PARAMETERS FOR A DUAL SCATTER TYPE LASER DOPPLER VELOCIMETER 
C ( LDV ) ♦ THE LDV MAY BE IN EITHER FORWARD SCATTER MODE OR BACK SCATTER 

C MODE* WHERE FORWARD SCATTER CONTAINS ANGLES FROM 0.0 DEGREES TO 90.0 
C DEGREES* AND BACK SCATTER CONTAINS ANGLES FROM 90.0 DEGREES TO 180.0 
C DEGREES. THE INPUT PARAMETERS NEEDED ARE* S I GNAL-TO-NO I SE RATIO AT 
C THE PHOTOMULT I PL IER* INPUT LASER POWER. FOCAL LENGTH OF INPUT 
C FOCUSING LENS. RECEIVER LENS DIAMETER. DISTANCE FROM SAMPLE VOLUME 
C TO RECEIVING LENS. LASER WAVELENGTH, LASER BEAM DIAMETER, STARTING 
C CROSS BEAM ANGLE. INCREMENT ANGLE, AND FINAL ANGLE, PARTICLE 
C DISTRIBUTION AND INDEX OF REFRACTION. 

C THE OUTPUT CONTAINS THE SIZE OF THE SAMPLE VOLUME, SCATTERED 

C LIGHT POWER RECEIVED, MIE INTENSITY FUNCTION, REQUIRED PARTICLE 
C DENSITY AND PHOTOTUBE OUTPUT CURRENT AND POWER. 

C SN = REQUIRED S I GNAL-TO-NO I SE RATIO IN DB . 

C PO = INPUT LASER POWER IN WATTS 
C FI = FOCAL LENGTH OF INPUT LENS IN METERS. 

C DR = DIAMETER OF RECEIVING LENS IN METERS. 

C Z = DISTANCE FROM SAMPlE VOLUME TO RECEIVING LENS IN METERS. 

C XL AM = LASER WAVELENGTH IN METERS. 

C D = INPUT LASER DIAMETFR IN METERS. 

C DTHETA = STARTING CROSF BEAM ANGLE IN DEGREES. 

C DETH = INCREMENT CROSS BEAM ANGLE IN DEGREES. 

C DTHELA = ENDING CROSS REAM ANGLE IN DEGREES. 

C S = STARTING PARTICLE DIAMETER IN MICRONS. 

C DELTA = INCREMENT PARTICLE DIAMETER IN MICRONS. 

C SLAST = ENDING PARTICLE DIAMETER IN MICRONS. 

C SREF = PARTICLE INDEX OF REFRACTION. 


READ 1. SN. PO. FI, DR , Z 
READ 2. XLAM. D 

READ 3, DTHETA, DETH, DTHELA 
READ 3. S, DELTA, SLAST 

READ 3, SREF 

1 FORMAT (F15.6) , 

2 FORMAT (E15.6) . 

3 FORMAT C3F10.8) 

SN9 = SN 

S9 = S 

DELTA9 = DELTA 
SLAST9 = SLAST 
P = PO/2.0 
1000 SN = SN9 
S = S9 

DELTA = DELTA9 
SLAST = SLAST9 
PRINT 100 
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100 FORMAT ( 1 HI »20X« *LDV DUAL SCATTER MODE, INPUT S/N*////) 


PRINT 

1 0 


PRINT 

1 1 , 

SN 

PRINT 

12, 

PO 

PRINT 

13, 

FI 

PRINT 

14, 

DTHFTA 

PRINT 

15, 

DR 

PR I NT 

16, 

Z 

PRINT 

17, 

XLAM 

PRINT 

18, 

D 

PRINT 

19, 

S, SLA 

PRINT 

20, 

SREF 


THETA = DTHETA*0. 01 745329252 

PIM = 3.I4159E-06 

s = PIM*S/XLAM 

DELTA = PIM*DELTA/XLAM 

SLAST - PIM*SLAST/XLAM 

H = 6.626E-34 , 

C = 3.0E08, 

DELF = 1 • 0E06 * 

B = 1.2, 

E = 1 . 6E — 1 9 
G = 8.0E+04 
XNU = 0.05 
ELI = 0.92, 

EL2 = 0*45, 

R = 50.0, 

SN = SN/10.0, 

SN * 10.0**SN, 

C CALCULATION OF REQUIRED SCATTERED POWER GOING TO PHOTOTUBE 

PS = 2. 0*H*C*DELF*B*gN/ ( XLAM*XNU ) 

PRINT 21 
PRINT 22, PS 

C CALCULATION OF REQUIRED SCATTERED POWER 

PSC = PS/ ( EL 1 *EL2 ) , 

PRINT 23, PSC 

C CALCULATION OF INPUT POWER DENSITY 

DA = 1 . 22*XLAM*F 1 /D 

PRINT 24, DA 

ACL = 3 • 1 4 1 59*DA*DA/4 .0 , 

PRINT 25, ACL 
F = P/ACL, 

PRINT 26, F 

C CALCULATION OF MIE INTENSITY FUNCTION 

CALL MIEABCS, DELTA, SLAST, SREF ) 

A A SUM = 0.0 
COUNT =0.0 

IF (DTHETA.GT.90.0) G O TO 201 
DTHETA = DTHETA/2.0 
GO TO 210 

201 DTHETA = 180.0 - DTHETA 

DTHETA = DTHETA/2.0 
DTHETA = 180.0 - DTHETA 

210 CONTINUE 

CALL M I E A C ( D THE T A , AASUM, COUNT) 

XIMIE = A A SUM /COUNT 
IF (DTHETA.GT.90.0) GO TO 202 
DTHETA = DTHETA*2.0 
GO TO 211 
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202 DTHETA = 180.0 - DTHfTA 

DTHETA = DTHET A*2 • 0 
DTHETA = 180.0 - DTHFTA 

211 CONTINUE 

PRINT 27, XIMIE 

C CALCULATION OF SAMPLE VOLUME , 

XL = 1 .41421 4*DA*C0S (THETA/2.0 )/SIN( THET A ) 

PRINT 34, XL. 

V = 3. 1 41 59*DA*DA *DA*COS( THET A/2.0 >*S I N ( THET A/2 . 0 ) / (3 . 0*S I N ( THETA ) 
1 *S I N (THET A ) ) 

PRINT 33, V 

C CALCULATION OF SOLID ANGLE , 

AR = 3. 141 59*DR*DR/4 .0, 

OMEGA = AR/(Z*Z) , 

PRINT 28, OMEGA 

C CALCULATION OF REQUIRED PARTICLE DENSITY 

XK = 2.0*3. 1 4159/XLAM, 

XNMIE = XK*XK*PSC/ ( F*V*0MEGA*X I M IF ) , 

XNMIE = XNMIE/10.0**6 , 

PRINT 29, XNMIE 

C CALCULATION OF OUTPUT POWER , 

PI = 2 • 0*PS 

XN = XLAM*P1 / < H*C > , 

PRINT 30, XN 

XIOUT = XN*G*E*XNU 

XIOUT = XIOUT*0.5857, 

PRINT 31 , XIOUT 

XPOUT =* X IOUT*X I 0UT*R* 1 000 .0, 

XPOUT = 1 0. 0*ALOG1 O ( XPOUT ) , 

PRINT 32, XPOUT 

IF (DTHETA. GE.DTHELA ) GO TO 200 
DTHETA = DTHETA + DFTH 
IF (DTHETA. GT. 90.0) GO TO 220 
THFTA2 = DTHETA*0 .0 1 745329252/2 . 0 
Z = DR/ (2 .0*TAN ( THETA2 ) ) 

GO TO 230 
220 Z = FI 
230 CONTINUE 

GO TO 1000 
200 CONTINUE 

10 FORMAT ( 20X , * INPUT PARAMETERS*//) 

11 FORMAT (* SIGNAL-TO-NOISE RATIO *F10.4* DB*/ ) 

12 FORMAT (* LASER POWER *F10.4* WATTS*/) 

13 FORMAT (* INPUT LENS FOCAL LENGTH *F10.4* METERS*/) 

14 FORMAT (* SCATTERING ANGLE *FI0.4* DEGREES*/) 

15 FORMAT (* COLLECTING LENS DIAMETER *F10.4* METERS*/) 

16 FORMAT (* DISTANCE FROM SAMPlE VOLUME TO COLLECTING LENS * 

1F10.4* METERS*/) 

17 FORMAT (* LASER WAVELENGTH *E15.6* METERS*/) 

18 FORMAT (* DIAMETER OF INPUT LASER BEAM *E15.6* METERS*/) 

19 FORMAT (* PARTICLE RANGE *F8.4* MICRONS TO*F8.4* MICRONS*/) 

20 FORMAT (* INDEX OF REFRACTION OF THE PARTICLES *F10.4* + J* 

1 F10.4////) 

21 FORMAT ( 20X , * OUTPUT PARAMETERS*//) 

22 FORMAT (* SCATTERED POWER RECEIVED BY PHOTOMULTIPLIER *E15.6 

2* WATTS*/) 

23 FORMAT (* SCATTERED POWER REQUIRED FROM FLOW *E15.6* WATTS*/) 

24 FORMAT (* AIRY DISK DIAMETER *E15.6* METERS*/) 

25 FORMAT (* CROSS SECTIONAL AREA OF SAMPLE VOLUME *E15.6 
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SAMPLE VOLUME *E15.6 


3* SQMETERS*/) 

26 FORMAT <* LASER POWER DENSITY AT 
4* WATTS/SQMETER*/ ) 

MIE INTENSITY FUNCTION *E15.6/> 

SOLID ANGLF OF THE VIEWING OPTICS *E15#6* SR*/) 

REQUIRED PARTICLE DENSITY *E15.6* PART I CLES/CC*/ ) 

INPUT POWER TO PHOTOMULTIPLIER *E15.6* PHOTONS/SEC*/ 

PHOTOMULT I PL I ER OUTPUT CURRENT *E15.6* AMPS*/) 

PHOTOMULTIPLIER OUTPUT POWER *F10*4* DBM*////) 

SAMPLE VOLUME *E15.6* CUBMETERS*/) 

LENGTH OF THE SAMPLE VOLUME *E15.6* METERS*/) 


27 

FORMAT 

( * 

28 

format 

(* 

29 

FORMAT 

(* 

30 

FORMAT 

(* 

5) 


31 

FORMAT 

(* 

32 

FORMAT 

(* 

33 

FORMAT 

(* 

34 

FORMAT 

(* 


STOP* 



END* 



SUBROUTINE M I E AB ( S * DELTA, SLAST * SREF ) 

C MIE A AND B MAIN PROGRAM, 

D I MENS ION SREF (2 ) i WN (2 ) * WNM1 (2 ) * WNM2 (2) » ALAST (2) * A NEXT ( 2 ) * REF (2 ) * A 
*(2 >* 8 < 2 ) * 

COMMON REF • X * A # 3 # FN « WN * WNM 1 « WNM2 tALAST » ANEXT * 

INTEGER ZERO 1 
ZERO 1 = 0 
LUM=9 

•1101 CONTINUE 

CALL RECOUT (LUN» 1*0# ZERO 1 # A ( 1 )*A(2)*3(1)*8(2)) 

CALL RECOUT ( LUN * 1 * 0 * F * SREF ( 1 ) * SREF ( 2 ) ) 

C INITIALIZATION* 

101 FN-O.O, 

T E ST = 1 • 2*S + 9 • * 

X = S , 

WN ( 1 )=SIN(X) i 
WN (2 ) =COS (X ) * 

WNM 1 ( 1 ) = WN (2 ) * 

WNM 1 (2 ) =-WN ( 1 ) , 

REF ( 1 ) = SREF ( 1 ) , 

REF ( 2 ) = SREF ( 2 ) , 

U^ REF (1 )*X* 

V = -RFF (2 >*X* 

SU = S IN (U ) * 

CU = COS CU ) , 

W-EXP ( V ) * 

Wl=l i/Wf 
SHV=0* 5* ( W-Wl ) * 

CHV-0.5* (W + Wl ) * 

T E M 1 = 1 ./ (SU*SU+SHV*SHV ) * 

ANFXT( 1 )=SU*CU*TEM1 , 

ANFXT ( 2 )=SHV*CHV*TEMI , 

C STEP UP FOR NEXT CALCULATION* 

3 ALAST ( 1 ) = ANEXT ( I ) * 

AL AST ( 2 ) = ANEXT ( 2 ) * 

WNM2 ( 1 ) — WNM 1 ( 1 ) * 

WNM2 (2 ) = WNM 1 (2 ) % 

WNM 1 ( 1 ) =WN ( 1 ) % 
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WNM1 (2 ) = WN ( 2 ) * 

FN=FN+ 1 • * 

C COMPUTE A AND B, 

CALL NEXT AB * 

N=FN • 

CALL RECOUT (LUN iltO<N<A(l ) , A C2 ) , B C 1 ) ,3 (2 > ) 

I F ( FN • GE • TEST ) GO TO 9* 

IF ( ( AB5 < A (1 ) ) • G T • 1 • E —8 ) • OR • ( ABS ( A (2 ) ) • GT • I .E-8 ) .OR. ( ABS <B ( 1 ) ) 

1 # GT • 1 . E— 8 ) • OR • ( ABS <R(2 ) ) .GT. 1 .E-8 ) ) GO TO 3 

9 IF<S.GE.SLAST)GO TO 1C. 

12 S = S + DELT A , 

GO TO 1101, 

10 CONTINUE 

700 CALL RECOUT (LUN ♦ 1 .0 .ZEROl , A ( 1 ) , A (2 ) .3 ( 1 ) *8 (2 ) ) 

END FILE 9 
901 REWIND 9 
RETURN 
END , 


SUBROUTINE NEXTAB 

D I MENS ION WN (2 ) , WNM 1 (2 ) , WNM2 (2) ♦ ALAST (2) .ANEXT (2 ) * REF (2),A(2),3C2) 
D I MENS I ON STORE 1(2). STORE 2 ( 2 ) ♦ ADENOM ( 2 ) * SDENOM C2 ) , ANUM (2 ) ♦ 3NUM ( 2 ) * 
COMMON REF ,XiA,B,FN, WN , WNM 1 , WNM 2 . AlAST .ANEXT , 

C CALCULATE NEXT A + B FROM |_AST, 

DENI = ( REF ( 1 ) *REF ( 1 ) + REF ( 2 ) *REF (2 ) }*Xi 
TERM1 = FN/DEN I * 

TERM2 = TERM 1 *REF (2 ) , 

TERM 1 =-TERMl *REF Cl), 

DENOM 1 = — TERM 1 -ALAST Cl). 

DENOM2 = — TERM 2 — AL A S T ( P ) , 

TFMP= 1 • / C DENOM 1 *DENOM 1 +DENOM2*DE N0M2 > • 

ANFXTC 1 )= TERM 1+DENOMl* TEMP, 

ANEXT (2 ) = TE R M2 — DENOM TEMP . 

FACT = ( 2 • *FN— 1 . ) /X . 

WN ( 1 )=FACT*WNM1 ( 1 )-WNM2 ( 1 ) . 

WNC2 )=FACT*WNM1 (2)-WNM2(2). 

STORE=FN/X. 

CALL DIVCPXC ANEXT. REF. ST0RE1 ). 

STORE l C 1 ) =STORE 1 ( 1 )+ STORE, 

CALL MPYCPX ( STORE 1 . WN . ST0RE2 ) . 

CALL SUBCPXCST0RE2.WNM1 .ADENOM) . 

ANUM ( 1 ) =WN ( 1 ) * ST ORE 1 C 1 ) . 

ANUM (2 ) =WN ( 1 ) *ST ORE 1 (Z ) . 

ANUM C 1 ) =ANUM ( 1 ) —WNM 1 Cl). 

CALL DIVCPXC ANUM. ADENOM, A). 

CALL MPYCPXCANEXT.REF.STOREl ) , 

ST ORE 1 ( 1 ) “ST ORE 1 ( 1 )+STORE, 

CALL MPYCPX ( STORE 1 , WN. ST0RE2 ) 

CALL SUBCPX ( STORE2, WmM 1 , 3DENOM ) , 

BNUM ( 1 ) = WN C 1 )*ST0RE1 Cl). 

BNUM (2 ) = WN ( 1 ) *ST ORE 1 (2 ) , 

BNUM C 1 ) =BNUM ( 1 ) -WNM1 C 1 ) . 
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CALL DIVCPX(BNUM.SDEMOM.B) , 

RETURN i 

END. 


SUBROUTINE Dl VCPX (X. Y.Z > . 
DIMENSION X(2 ) . Y (2 ) .7 ( 2 > ♦ 

DENOM = Y( 1 )*Y ( 1 )+Y (2 )*Y<2 > , 

Zi = < X ( 1 >*Y ( 1 ) + X (2 >*Y (2 > ) /DENOM . 
Z2 = C X ( 2 ) * Y ( 1 )-X ( 1 ) * Y (2 ) ) /DENOM. 
Z( 1 )=Z1 . 

Z (2 ) =Z2 . 

RETURN. 

END. 


SUBROUT I NE MPYCPX (X. Y.Z ) . 
D I MENS ION X ( 2 ) . Y ( 2 ) . 7 ( 2 ) . 
Zl =X ( 1 ) * Y ( 1 )-X (2 )*Y (? ) . 

Z( 1 )=Z1 . 

Z ( 2 )=Z2. 

RETURN. 

END, 


SUBROUTINE SUBCPX(X. Y.Z ) . 
DIMENSION X(2).Y(2).7(2). 
Z ( I ) =X ( 1 )-Y ( 1 ) , 

Z ( 2 ) =X (2 )-Y ( ? ) , 

RETURN, 

FND, 


SUBROUTINE M I E AC (DTHETA, AASUM. COUNT) 

** PROGRAM TO CALC. QE.QS, AND 1(180) USING DE I RMENDJ I AN— S A-S AND B~S 
C 

DIMENSION A ( 500 .2 ) , B (500 . 2 ) , M < 500 ) , P I < 500 ) 4 T AU ( 500 ) 

C 

1000 1=1 

LUN = 9 

THETA=DTHETA*.0 1 7453^92519943 
STHET A = S I N ( THETA ) 
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STHESQ=STHETA**2 
CTHETA=COS( THETA ) 

CTHESQ=CTHETA**2 

PI (1 ) = 1 • 0 

PI (2 ) = 3 • *C THET A 

TAU( 1 ) =CTHET A 

T AU (2 ) =3 • * ( CTHESQ-STHESQ > 

CALL REC IN (LUN « 1 • 0* M ( I ).A(I#1)*A(I,2),B(I«1 ) «B( I ,2 ) ) 

5 CALL RECIN (LUN ♦ 1 ♦ 0 ♦ ALPHA ♦ ET AR * ETA I ) 

I F (EOF « 9 ) 90 1 « 1 

1 DO 2 11 = 1 * 250 

r = i r 

CALL REC IN (LUN *l*0.M<I).A(l.l)*A(I,2)«B<I%l)tB(It2>) 

IF (M ( I ) ) 2 . 3 ♦ 2 

2 CONTINUE 

3 CONTINUE 

MAX =1-1 
SUM 1 = 0.0 

SUM2 = 0.0 
C=-l .0 

DO 11 1=1 * MAX 

X I = I 

SUM 1 = (2.*Xl+l .0 >* ( A( I « 1 )+B ( I ♦ 1 ) ) + SUM 1 

SUM2 = ( 2 • *X I + 1 • 0 ) * ( A ( I < 1 )**2+A (1,2 )**2+B ( l , 1 ) **2 + B (1,2 )**2 )+SUM2 

11 CONTINUE 

QF = 2.0*SUM1/(ALPHA**2) 

QS = 2. 0*SUM2/ ( ALPHA**2 ) 

IF (DTHETA, NE •180# > GO TO 12 

SUM4=0 # 0 

SUMRE = 0*0 

SUMIM = 0,0 

DO 4 1=1 ♦ MAX 

C=C* (-1 .0 ) 

XI = I 

SUMRE = (Xl+0.5)*(A<I.l) — B (1*1)) *C+ SUMRE 

SUMIM = (Xl+0.5)*(A( 1*2)— B( I ) ) *C+SUMI M 

4 CONTINUE 

SOMGRE = SUMRE*SUMRF 
SQMGIM = SUM I M*SUMI M 
SUM3 = SQMGRE 4- SQMGIM 
GO TO 30 

12 CONTINUE 

DO 13 NN = 3 t MAX 
XN = NN 

TWONM1 =2 • *XN— 1 . 

PI (NN ) = (TW0NM1 *PI ( NN — 1 ) *CTHETA-XN*P I (NN-2 ) )/(XN- 1 • ) 

TAU(NN )=CTHETA* (PI (NN)~PI (NN-2) )-(TWONM! >*STHESQ*PI (NN-1 >+TAU(NN-2 
1 ) 

13 CONTINUE 
SS SRI =0.0 
SSSI 1=0.0 
SSSR2=0 • 0 
SSS I 2 = 0 • 0 

DO 14 NN=liMAX 
XN = NN 

TEMP =(2.*XN+1 • )/(XN* ( XN+ 1 . ) ) 

TEMPP I =TEMP*P I (NN ) 

TEMPT AU= TEMP* T AU (NN ) 

SSSR1 =TEMPP I *A (NN# 1 ) +TEMPT AU*B ( NN , 1 J+SSSR1 
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SSSI 1 =TEMPPI*A (NN*2 ) +TEMPT AU*B { NN • 2 )+SSS I 1 
SSSR2-TEMPTAU*A (NN* 1 ) + TEMPP I *B ( NN ♦ 1 >+SSSR2 
SSSI2=TEMPTAU*A ( NN * 2 ) +TEMPP 1 *B ( NN * 2 )+SSS 12 
14 CONTINUE 

5UM3 = SSSR1**2 + SSS I 1**2 
SUM4=SSSR2**2+SSS 12**2 
AASUM = SUM4 + AASUM 
COUNT = COUNT + 1 *0 
30 CONTINUE* 

GO TO 5 
C 

901 CONTINUE* 

700 CONTINUE 
REWIND 9 
RETURN 
END 

-00 


15.0 
0.05 
0.6096 
0.0762 
4 *366 

6328000E-07 
2000000E— 03 

170.0 0.25 179.0 

0.4 0.005 2.0 

1 .48 
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BEAM SPLITTER 




CD 



Figure 2.- Schematic of cross-beam laser Doppler velocimeter in reference-scatter mode. 



PHOTOMULTIPLIER POSITION 
FOR BACK-SCATTER MODE 


Figure 3.- Schematic of cross-beam laser Doppler velocimeter in dual-scatter — forward-scatter 

and dual-scatter — back-scatter modes. 
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Figure 7.- Comparison of signal-to-noise ratio determined from a seeding 
density of 1000 partieles/cm 3 as a function of cross-beam angle for the 
three LDV modes using the 5145 A Argon laser line at 1.0 W and a 
noise bandwidth of 10® Hz. 



OUTPUT CURRENT , 



ANGLE BETWEEN CROSS BEAMS , deg 

Figure 8.- Comparison of photomultiplier output current determined from a seeding 
density of 1000 particles/cm^ as a function of cross-beam angle for the three 
LDV modes using the 5145 A Argon laser line at 1.0 W. 




Figure 9.- Comparison of photomultiplier output power determined from a seeding 
density of 1000 particles/cm^ as a function of cross-beam angle for the three 
LDV modes using the 5145 A Argon laser line at 1.0 W. 
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PARTICLE DENSITY, particles/c 



ANGLE BETWEEN CROSS BEAMS, deg 

Figure 10,- Determination of required seeding density as a function of cross-beam 
angle for the three LDV modes using a signal-to-noise ratio of 15 dB, a noise 
bandwidth of 10 6 Hz, and 5145 A Argon laser line at 1.0 W. 



PARTICLE DENSITY , particles/. 



Figure 11.- Determination of required seeding density as a function of cross-beam 
angle for the dual-scatter — forward-scatter mode and the reference-scatter 
mode using a 50-mW He/Ne laser and a 1.0-W Argon laser at 5145 A for a 
signal -to-noise ratio of 15 dB and a noise bandwidth of 10® Hz. 
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INPUT LASER POWER, W 


Figure 12.- Calculation of input laser power required to obtain molecular scatter 
needed for a signal-to-noise ratio of 15 dB using a 5145 A Argon laser line as 
a function of molecular density, with a noise bandwidth of 10^ Hz and an S-20 
photomultiplier response. 
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MOLECULAR DENSITY, molecules/cm' 



INPUT LASER POWER, W 


Figure 14.- Calculation of required input laser power to obtain molecular scatter 
needed for a signal-to-noise ratio of 15 dB as a function of molecular density 
for several laser wavelengths in air, with a noise bandwidth of 10^ Hz and an 
S-20 photomultiplier response. 
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